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If the scale of quantum gravity is near a TeV, the LHC will be producing one black hole (BH) 
about every second. The BH decays into prompt, hard photons and charged leptons is a clean 
signature with low background. The absence of significant missing energy allows the reconstruction 
of the mass of the decaying BH. The correlation between the BH mass and its temperature, deduced 
from the energy spectrum of the decay products, can test experimentally the higher dimensional 
Hawking evaporation law. It can also determine the number of large new dimensions and the scale 
of quantum gravity. 
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Introduction: An exciting consequence of TeV-scale 
quantum gravity M is the possibility of production of 
black holes (BHs) |§ |, Q at the LHC and beyond. The 
objective of this paper is to point out the experimental 
signatures of BH production. Black holes are well under- 
stood general-relativistic objects when their mass Mbh 
far exceeds the fundamental (higher dimensional) Planck 
mass Mp ~TeV. As Mbh approaches Mp, the BHs be- 
come "stringy" and their properties complex. This raises 
an obstacle to calculating the production and decay of 
light BHs, those most directly accessible to the LHC, 
where the center-of-mass (c.o.m.) energy of colliding 
beams is comparable to the Planck mass. In what follows, 
we will ignore this obstacle and estimate the properties 
of light BHs by simple semiclassical arguments, strictly 
valid for Mbh 3> Mp. We expect that this will be an 
adequate approximation, since the important experimen- 
tal signatures rely on two simple qualitative properties: 
(i) the absence of small couplings and (ii) the "demo- 
cratic" (flavor independent) nature of BH decays, both 
of which may survive as average properties of the light 
descendants of black holes. Nevertheless, because of the 
unknown stringy corrections, our results are approximate 
estimates. For this reason, we will not attempt selective 
partial improvements - such as time dependence, angular 
momentum, charge, hair, and other higher-order general 
rclativistic refinements - which, for light BHs, may be 
masked by larger unknown stringy effects. We will focus 
on the production and sudden decay of Schwarzschild 
black holes. 

Production: The Schwarzschild radius Rs of an (4 + 
n)-dimensional black hole is given by || : 
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assuming that extra dimensions are large (3> Rs)- 

Consider two partons with the c.o.m. energy Vs = 
Mbh moving in opposite directions. Semiclassical rea- 
soning suggests that if the impact parameter is less than 
the (higher dimensional) Schwarzschild radius, a BH with 
the mass Mbh forms. Therefore the total cross section 
can be estimated from geometrical arguments H, and is 
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(see Fig. ||a). 

This expression contains no small coupling constants; 
if the parton c.o.m. energy y/l reaches the fundamental 
Planck scale Mp ~ TeV then the cross section if of or- 
der TeV~ 2 « 400 pb. At the LHC, with the total c.o.m. 
energy y/s = 14 TeV, BHs will be produced copiously. 
To calculate total production cross section, we need to 
take into account that only a fraction of the total c.o.m. 
energy in a pp collision is achieved in a parton-parton 
scattering. We compute the full particle level cross sec- 
tion using the parton luminosity approach (after Ref. M): 

da{pp -> BH + X) dL 
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where the parton luminosity dL / dM B n is defined as the 
sum over all the initial parton types: 
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and fi(xi) are the parton distribution functions (PDFs). 
We used the MRSD— ' [§ PDF set with the Q 2 scale 
taken to be equal to Mbh, which is within the allowed 
range for this PDF set, up to the LHC kinematic limit. 
The dependence of the cross section on the choice of PDF 
is ~ 10%, i.e. satisfactory for the purpose of this esti- 
mate. 

The differential cross section dtj/dMsH for the BH 
produced at the LHC is shown in Fig. |l|b for several 
choices of Mp. The total production cross section at 
the LHC for BH masses above Mp ranges from 0.5 nb 
for M P = 2 TeV, n = 7 to 120 fb for M P = 6 TeV and 
n = 3. If the fundamental Planck scale is ~ 1 TeV, LHC, 
with the peak luminosity of 30 fb _1 /year will produce 
over 10 7 black holes per year. This is comparable to the 
total number of Z's produced at LEP, and suggests that 
we may do high precision studies of TeV BH physics, as 
long as the backgrounds are kept small. 

Decay: The decay of the BH is governed by its Hawk- 
ing temperature T#, which is proportional to the inverse 
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FIG. 1: a) Parton- level production cross section, b) differen- 
tial cross section da/dMsn at the LHC, c) Hawking temper- 
ature, and d) average decay multiplicity for a Schwarzschild 
black hole. The number of extra spatial dimensions n = 4 is 
used for a)-c). The dependence of the cross section and Hawk- 
ing temperature on n is weak and would be hardly noticable 
on the logarithmic scale. 



radius, and given by [0: 



T H =M P 



M P 




(3) 



(see Fig. |l]b). As the parton collision energy increases, 
the resulting black hole gets heavier and its decay prod- 
ucts get colder. 

Note that the wavelength A = ^ corresponding to 
the Hawking temperature is larger than the size of the 
black hole. Therefore, the BH is, to first approxima- 
tion, a point-radiator and therefore emits mostly s- waves. 
This indicates that it decays equally to a particle on the 
brane and in the bulk, since it is only sensitive to the 
radial coordinate and does not make use of the extra an- 
gular modes available in the bulk. Since there are many 
more particles on our brane than in the bulk, this has the 
crucial consequence that the black hole decays visibly to 
standard model (SM) particles [Q, ||. 

The average multiplicity of particles produced in the 
process of BH evaporation is given by: (N) = / \ , 
where E is the energy spectrum of the decay products. 
In order to find (N), we note that the BH evaporation 
is a blackbody radiation process, with the energy flux 
per unit of time given by Planck's formula: 4^ ~ e *±c > 
where x = E/Th, and c is a constant, which depends on 
the quantum statistics of the decay products (c — — 1 for 
bosons, +1 for fermions, and for Boltzmann statistics). 



The spectrum of the BH decay products in the massless 
particle approximation is given by: 4^ ~ TiIie ~ J°±c ' 
In order to calculate the average multiplicity of the par- 
ticles produced in the BH decay, we use the average of 
the distribution in the inverse particle energy: 
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where a is a dimensionless constant that depends on the 
type of produced particles and numerically equals 0.68 for 
bosons, 0.46 for fermions, and ^ for Boltzmann statistics. 
Since a mixture of fermions and bosons is produced in 
the BH decay, we can approximate the average by using 
Boltzmann statistics, which gives the following formula 



for the average multiplicity: (N) 
for Hawking temperature, we obtain: 
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Using Eq. 
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Eq. (|J) is reliable when the mass of the BH is much 
larger than the Hawking temperature, i.e. (N) ^> 1; oth- 
erwise, the Planck spectrum is truncated at£w Mbh/2 
by the decay kinematics 10 . The average number of 



particles produced in the process of BH evaporation is 
shown in Fig. [j]d, as a function of M^n/Mp, for several 
values of n. 

We emphasize that, throughout this paper, we ignore 
time evolution: as the BH decays, it gets lighter and hot- 
ter and its decay accelerates. We adopt the "sudden ap- 
proximation" in which the BH decays, at its original tem- 
perature, into its decay products. This approximation 
should be reliable as the BH spends most of its time near 
its original mass and temperature, because that is when 
it evolves the slowest; furthermore, that is also when it 
emits the most particles. Later, when we test the Hawk- 
ing mass-temperature relation by reconstructing Wien's 
dispacement law, we will minimize the sensitivity to the 
late and hot stages of the BHs life by looking at only 
the soft part of the decay spectrum. Proper treatment 
of time evolution, for Mbh ~ Mp, is difficult, since it 
immediately takes us to the stringy regime. 

Branching Fractions: The decay of a BH is ther- 
mal: it obeys all local conservation laws, but otherwise 
does not discriminate between particle species (of the 
same mass and spin). Theories with quantum gravity 
near a TeV must have additional symmetries, beyond 
the standard SU(3) x SU(2) x U(\), to guarantee pro- 
ton longevity, approximate lepton number (s) and flavor 
conservation ]TH |. There are many possibilities: discrete 
or continuous symmetries, four dimensional or higher di- 
mensional "bulk" symmetries jl2j . Each of these possi- 
ble symmetries constrains the decays of the black holes. 
Since the typical decay involves a large number of par- 
ticles, we will ignore the constraints imposed by the few 
conservation laws and assume that the BH decays with 
roughly equal probability to all off « 60 particles of the 
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SM. Since there are six charged leptons and one photon, 
we expect ~ 10% of the particles to be hard, primary lep- 
tons and ~ 2% of the particles to be hard photons, each 
carrying hundreds of GeV of energy. This is a very clean 
signal, with negligible background, as the production of 
SM leptons or photons in high-multiplicity events at the 
LHC occurs at a much smaller rate than the BH produc- 
tion (see Fig. ||). These events are also easy to trigger on, 
since they contain at least one prompt lepton or photon 
with the energy above 100 GeV, as well as energetic jets. 

Test of the Hawking's radiation: Furthermore, 
since there are three neutrinos, we expect only ~ 5% 
average missing transverse energy ($ T ) per event, which 
allows us to precisely estimate the BH mass from the 
visible decay products. We can also reconstruct the BH 
temperature by fitting the energy spectrum of the decay 
products to the Planck's formula. Simultaneous knowl- 
edge of the BH mass and its temperature allows for a test 
of the Hawking's radiation and can provide an evidence 
that the observed events come from the production of 
BH, and not from some other new physics. 

There are a few important experimental techniques 
that we will use to carry out the numerical test. First 
of all, to improve precision of the BH mass reconstruc- 
tion we will use only the events with $ T consistent with 
zero. Given the small probability for a BH to emit a 
neutrino or a graviton, total statistics won't suffer ap- 
preciably from this requirement. Since BH decays have 
large jet activity, the Mbh resolution will be dominated 
by the jet energy resolution and the initial state radiation 
effects, and is expected to be ~ 100 GeV for a massive 
BH. Second, we will use only photons and electrons in the 
final state to reconstruct the Hawking temperature. The 
reason is twofold: final states with energetic electrons 
and photons have very low background at high Vs, and 
the energy resolution for electrons and photons remains 
excellent even at the highest energies achieved in the pro- 
cess of BH evaparation. We do not use muons, as their 
momenta are determined by the track curvature in the 
magnetic field, and thus the resolution deteriorates fast 
with the muon momentum growth. We also ignore the 
r-lepton decay modes, as the final states with r's have 
much higher background than inclusive electron or pho- 
ton final states, and also because their energies can not 
be reconstructed as well as those for the electromagnetic 
objects. Fraction of electrons and photons among the fi- 
nal state particles is only ~ 5%, but the vast amount of 
BHs produced at the LHC allows us to sacrifice the rest of 
the statistics to allow for a high-precision measurement. 
(Also, the large number of decay particles enhances the 
probability to have a photon or an electron in the event.) 
Finally, if the energy of a decay particle approaches the 
kinematic limit for pair production, Mbh/2, the shape 
of the energy spectrum depends on the details of the BH 
decay model. In order to eliminate this unwanted model 
dependence, we use only the low part of the energy spec- 
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FIG. 2: Number of BHs produced at the LHC in the electron 
or photon decay channels, with 100 fb -1 of integrated lumi- 
nosity, as a function of the BH mass. The shaded regions 
correspond to the variation in the number of events for n be- 
tween 2 and 7. The dashed line shows total SM background 
(from inclusive Z(ee) and direct photon production). The 
dotted line corresponds to the Z(ee) + X background alone. 

trum with E < Mbh/2. 

The experimental procedure is straightforward: we se- 
lect the BH sample by requiring events with high mass (> 
1 TeV) and mutiplicity of the final state (N > 4) , which 
contain electrons or photons with energy > 100 GeV. We 
smear the energies of the decay products with the reso- 
lutions typical of the LHC detectors. We bin the events 
in the invariant mass with the bin size (500 GeV) much 
wider than the mass resolution. The mass spectrum of 
the BHs produced at the LHC with 100 fkr 1 of inte- 
grated luminosity is shown in Fig. [2] for several values of 
Mp and n. Backgrounds from the SM Z(ee)+ jets and 
7+ jets production, as estimated with PYTHIA JI3], are 
small (sec figure). 

To determine the Hawking temperature in each A/bh 
bin, we perform a maximum likelihood fit of the energy 
spectrum of electrons and photons in the BH events to 
the Planck formula (with the coefficient c determined by 
the particle spin), below the kinematic cutoff (Mbh/2). 
We then use the measured Mbh vs. Tjj dependence 
and Eq. (||) to determine the fundamental Planck scale 
Mp and the dimensionality of space n. Note that to 
determine n we can also take the logarithm of both sides 
of Eq. (§): 

log(T H ) = — l — log(M BH ) + const, (6) 
n + 1 

where the constant does not depend on the BH mass, but 
only on Mp and on detailed properties of the bulk space, 
such as shape of extra dimensions. Therefore, the slope of 
a straight-line fit to the log(T# ) vs. log(MBH) data offers 
a direct way of determining the dimensionality of space. 
This is a multidimensional analog of Wien's displacement 
law. Note that Eq. (J6|) is fundamentally different from 
other ways of determining the dimensionality of space- 
time, e.g. by studying a monojet signature or a virtual 
graviton exchange processes, also predicted by theories 
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FIG. 3: Determination of the dimensionality of space via 
Wien's displacement law at the LHC with 100 fb" 1 of data. 

TABLE I: Determination of Mp and n from Hawking's radi- 
ation. The two numbers in each column correspond to frac- 
tional uncertainty in Mp and absolute uncertainty in n, re- 
spectively. 

M P 1 TeV 2 TeV 3 TeV 4 TeV 5 TeV 
n = 2 1%/0.01 l%/0.02 3.3%/0.10 16%/0.35 40%/0.46 
n = 3 1%/Q.01 1.4%/0.06 7.5%/0.22 30%/1.0 48%/1.2 
n = 4 1%/0.01 2.3%/0.13 9.5%/0.34 35%/1.5 54%/2.0 
n = 5 l%/0.02 3.2%/0.23 17%/1.1 
n = 6 l%/0.03 4.2%/0.34 23%/2.5 Fit fails 
n = 7 l%/0.07 4.5%/0.40 24%/3.8 



mination of the mass and the temperature of the black 
hole, and may lead to a test of Hawking's radiation. 

It is desirable to improve our primitive estimates, es- 
pecially for the light black holes (Mbh ~ Mp); this will 
involve string theory. Nevertheless, the most telling sig- 
natures of BH production - large and growing cross sec- 
tions; hard leptons, photons, and jets - emerge from qual- 
itative features that are expected to be reliably estimated 
from the semiclassical arguments of this paper. 

Perhaps black holes will be the first signal of TeV-scale 
quantum gravity. This depends on, among other factors, 
the relative magnitude of Mp and the (smaller) string 
scale Ms- For Ms <C Mp, the vibrational modes of the 
string may be the first indication of the new physics. 

Note added: After the completion of this work, a re- 
lated paper jl5| has appeared in the LANL archives. 

Acknowledgments: We would like to thank Gia 
Dvali, Veronika Hubeny, Nemanja Kaloper, Elias Kir- 
itsis, Konstantin Matchev, and Lenny Susskind for valu- 
able conversations, and many of the participants of the 
"Physics at TeV Colliders" and "Avatars of M-Theory" 
workshops for their interest. 



with large extra dimensions. 

Test of the Wien's law at the LHC would provide a 
confirmation that the observed e + X and 7 + X event 
excess is due to the BH production. It would also be 
the first experimental test of the Hawking's radiation hy- 
pothesis. Figure || shows typical fits to the simulated 
BH data at the LHC, corresponding to 100 fb _1 of in- 
tegrated luminosity, for the highest fundamental Planck 
scales that still allow for determination of the dimension- 
ality of space with reasonable precision. The reach of the 
LHC for the fundamental Planck scale and the number 
of extra dimensions via Hawking's radiation extends to 
Mp ~ 5 TeV and is summarized in Table jlj] . 

Note, that the BH discovery potential at the LHC is 
maximized in the e/fi+X channels, where background is 
much smaller than that in the 7 + AT channel (see Fig. |^). 
The reach of a simple counting experiment extends up to 
Mp » 9 TeV (n = 2-7), where one would expect to see 
a handful of BH events with negligible background. 

Summary: Black hole production at the LHC may be 
one of the early signatures of TeV-scale quantum gravity. 
It has three advantages: 

Large Cross Section. Because no small dimension- 
less coupling constants, analogous to a, suppress the pro- 
duction of BHs. This leads to enormous rates. 

Hard, Prompt, Charged Leptons and Photons. 
Because thermal decays are flavor-blind. This signature 
has practically vanishing SM background. 

Little Missing Energy. This facilitates the deter- 
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Results presented at the Les Houches Work- 
shop "Physics at the TeV Colliders" (May 30, 
2001) and the "Avatars of M- Theory" confer- 
ence, ITP at Santa Barbara (June 7, 2001 ). 



http: / /online, itp. ucsb.edu/onlinc/mthcory_c01/ dimopoulos 

[1] N. Arkani-Hamed, S. Dimopoulos, and G. Dvali, Phys. 
Lett. B 429, 263 (1998); I. Antoniadis, N. Arkani-Hamed, 
S. Dimopoulos, and G. Dvali, Phys. Lett. B 436, 257 
(1998); N. Arkani-Hamed, S. Dimopoulos, and G. Dvali, 
Phys. Rev. D 59, 086004 (1999). 

[2] P.C. Argyres, S. Dimopoulos, and J. March-Russell, 
Phys. Lett. B441, 96 (1998V 

[3] T. Banks and W. Fischler, ^ep-th/990603£ . 

[4] R. Emparan, G.T. H orowitz, and R.C . Myers, Phys. Rev. 
Lett. 85, 499 (2000) |hep-th/000311£ ]. 

[5] R.C. Myers and M.J. Perry, Ann. Phys. 172, 304 (1986). 

[6] In fact the cross section is somewhat enhanced by initial- 
state attraction 4L 



[7] V.D. Barger and R.J.N. Phillips, "Collider Physics, Up- 
dated Edition," Addison- Wesley, 1997, p. 157. 

[8] A.D. Martin, R.G. Roberts, and W.J. Stirling, preprint 
RAL-92-078 (1992). 

[9] L. Susskind, private communication. 
[10] To avoid this limitation when performing numerical tests 
of Wien's law, we truncate the integrals in Eq. (Q) at the 
kinematic limit (Mbh/2). 
[11] Any form of new physics at a TeV requires such new 
symmetries. An example, for TeV-scale supersymmetry, 
is R-parity. 

[12] N. Arkani-Hamed and S. Dimopoulos, |iep-ph/9811353 . 

[13] T. Sjostrand, Comp. Phys. Comm. 82, 74 (1994). 

[14] We checked the effect of introducing an additional cut- 
off at E < Ms, where string dynamics is expected 
to affect the shape of the Planck spectrum, by taking 
Ms = Mp/2. As a result, the uncertainties in Mp and n 
increased by about a factor of 2. 

[15] S.B. Giddings and S. Thomas, |hep-ph/01062ig| , v2. 



